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Abstract 

We give a variant of the homogeneous Buchberger algorithm for positively graded lattice ide- 
als. Using this algorithm we solve the Sullivant computational commutative algebra challenge 1 . 

1 Introduction 

Suppose that / is a homogeneous ideal in a polynomial ring R = k[xo, . . . ,x n ] over a field k. The 
usual homogeneous Buchberger algorithm builds a Grobner basis for / by successively constructing 
truncated Grobner bases of increasing degrees. Suppose that / is saturated i.e. I = I = {geR\ 
(xq ■ ■■x n ) m g G 7,m ^> 0}. If we encounter a polynomial / divisible by a variable in degree d of 
the homogeneous Buchberger algorithm, then we may conclude that / reduces to zero modulo the 
already constructed truncated Grobner basis in degree < d for /. This simple observation also allows 
for detection of non-saturated ideals in some cases. 

Sullivant's challenge is about deciding if a specified set B of 145,512 binomials generate the 
kernel P of the (toric) ring homomorphism (p : k[xijk\ — > k[u{j, v^, Wjk] given by 

(p{x ijk ) UijVikWjk- 

where 1 < z, j, k < 4. We give a version of the homogeneous Buchberger algorithm with a Gebauer- 
Moller criterion specifically tailored to positively graded lattice ideals. Using an implementation 
of this algorithm in the software package GLATWALK 2 we deduce that the ideal J generated by B is 
strictly contained in P by showing that J cannot be saturated. In fact, we exhibit a specific binomial 
b of degree 14 in J\ J. 

1 http : //math . berkeley . edu/~seths/ ccachallenge . html 
7 http : //home . imf . au . dk/niels/GLATWALK 
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I am grateful to B. Sturmfels for stimulating my interest in Sullivant's computational commuta- 
tive algebra challenge. R. Hemmecke has made me aware that he and P. Malkin already computed 
the full Grobner basis of J using new algorithms in a new version of 4ti2 thereby answering Sulli- 
vant's challenge. In fact they prove that the "missing" binomials in Sullivant's challenge have degree 
14 and form an orbit under the action of a certain symmetry group. I am grateful to Hemmecke for 
verifying that b lies in this orbit. 

2 Preliminaries 

We let R = k[x\, ...,x n ] denote the ring of polynomials over a field k. We assign degrees to the 
variables by deg(xi) = a\, . . . ,deg(x n ) = a n , where a\,...,a n are positive integers. A monomial 

x v E R has degree deg(^: v ) = v\a\ H hv„a„, where v = (vi,...,v„). This gives the (positive) 

grading 

R = ® s >oR s , 

where R s = span k {x v | deg(x v ) = s}. For a monomial order -< on R and a subset S C R we let 
in^(S) = {in^(/) | / E S\ {0}}. A Grobner basis for an ideal I C R is a finite subset G El, such 
that (in x (G)) =in^(7). 

2.1 Truncated Grobner bases 

For a homogeneous ideal I in R and d E N we let 

1<A = ®Is- 

A J-truncated Grobner basis for / is a finite subset G <( i C /<^, such that (in^(G<^))<rf = in(7) <£ / 
i.e. we require only match of initial ideals up to degree d. Using the division algorithm it is easy to 
show that / G I <c i reduces to zero modulo the polynomials in a J-truncated Grobner basis for I. 

3 The homogeneous Buchberger algorithm with sat-reduction 

We call an ideal I saturated if I = I = {g e R \ (xq ■ ■ -x n ) m g E 7,m 3> 0}. This means that mf E I 
implies f El, where m is a monomial and / a polynomial in R. Let -< be a term order on R. For a 
polynomial f E R we let sat(/) denote / divided by the greatest common divisor of the monomials 
in /. We say that / sat-reduces to h modulo g if either h = sat(/) and deg(/z) < deg(/) or / reduces 
to h modulo g in the usual sense i.e. in^ (g) divides a term t in f and 

h = f-{t/m < {g))g. 

Notice that if / sat-reduces to h modulo g and f,g belongs to a saturated ideal I, then h El. A 
remainder in the division algorithm of / modulo a set of polynomials G using sat-reduction in each 
step is denoted / G ( sat ) . 

The 5-polynomial of two homogeneous polynomials is homogeneous of degree no less than 
the degrees of the polynomials. The (usual) reduction of a homogeneous polynomial of degree d 
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modulo a set of homogeneous polynomials gives a homogeneous polynomial of degree d. These 
observations give the homogeneous Buchberger algorithm as explained in (H), Theorem 11). We 
tailor the homogeneous Buchberger algorithm to the special case where input consists of a set B = 
{/i, . . . ,f r } C R of homogeneous polynomials generating a saturated ideal. This has the consequence 
that reduction of a homogenous polynomial / of degree d divisible by a variable is not necessary, 
since f/xj G I <t j reduces to zero using the already computed J-truncated Grobner basis G <c j. 

Algorithm 3.1 (Homogeneous Buchberger algorithm for saturated ideals) 

INPUT: Term order -<. Homogeneous polynomials B = ,f r } C R generating a saturated 

homogeneous ideal /. 

OUTPUT: Homogeneous polynomials G = {g\, . . .,g s } such that {g\, . . .,g s } is a minimal Grobner 
basis over -< for the ideal generated by B. 

(i) S pairs :=0; G:=0; 

(ii) while (B ^ or Spairs ^ 0) do 

(a) Extract 3 a polynomial / of minimal degree in B U Spairs. 

(b) Compute g := y G ( sat ) 5 continue if the degree drops in a sat-reduction step in the division 
algorithm; 

(c) if (g = 0) continue; 

(d) G:=GU{g}; 

(e) Append 5-polynomials S(g, h) to Spairs for every hEG\ {g}. 
Remark 3.2 

(i) After step (liidt in Algorithm 13.11 the polynomials of degree < d in G form a minimal d- 
truncated Grobner basis of /, where d is the minimal degree of the polynomials in B U Spairs. 

An easy modification to algorithm d3.ll) may detect if / is not saturated. If the sat-reduction 

yG(sat) of y 

is non-zero and has lower degree than /, then we may deduce the existence of a monomial 
x v and a polynomial g such that x v g G /, but g^Ll. 

Algorithm 3.3 (Homogeneous Buchberger algorithm with saturation check) 
INPUT: Homogeneous polynomials B — {/i, . . . ,f r } and a term order -<. 

OUTPUT: Homogeneous polynomials G = {gi, . . .,g s } such that {gi, . . .,g s } is a minimal Grobner 
basis over -< for the ideal / generated by B or proof that / is not saturated. 

(i) Spairs :=0; G:=0; 
3 This means that / is deleted from the relevant list after it is extracted 
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(ii) while (B ^ or Spairs ^ 0) do 

(a) Extract a polynomial / of minimal degree d in B U Spairs. 

(b) g:=/ G ( sat ); 

(c) if (g = 0) continue; 

(d) if(deg(g) < d) 

(i) OUTPUT / as proof that / is not saturated and HALT 

(e) G:=GU{g}; 

(f) Append 5-polynomials S(g, h) to Spairs for every h G G\ {g}. 
Example 3.4 We give a simple example illustrating algorithm (13 .3b . 

(i) Consider the input B = {xz — y 2 ,x 4 — y 3 } along with the reverse lexicographic term order x -< 

v -<; z. 

(ii) The ideal / generated by B is homogeneous in the grading deg(*) = 3,deg(y) = 4,deg(z) = 5 
and deg(*z — y 2 ) = 8 < deg(x 4 — y 3 ) = 12. 

(iii) After the first loop we have B = {y 3 — x 4 }, G = {y 2 — xz} and Spairs = 0, where G is a 12- 
truncated Grobner basis of /. 

(iv) In the second loop we sat-reduce y 3 —x 4 modulo y 2 — xz and get yz — x 3 . As deg(yz — x 3 ) = 
9<deg(y 3 -x 4 ) = 12, we conclude that / is not saturated. 

(i) Now suppose that B = {y 2 — xz, yz — x 3 } in the same grading. 

(ii) After the second loop we have 

5 = 

G = {y 2 -xz,yz-x 3 } 
Spairs = {yx 3 — z 2 x}, 
where G is a 13-truncated Grobner basis of /. 

(iii) Now yx 3 — z 2 x sat-reduces to z 2 — yx 2 modulo G. We conclude that / is not saturated. 

(i) Now proceed with B = {y 2 — xz, yz — x 3 ,z 2 — yx 2 } . 

(ii) After a few loops we have 

5 = 

G = {y 2 -xz,yz-x 3 ,z 2 -yx 2 } 
Spairs = {y 2 x 2 — zx 3 }, 

where G is a 14-truncated Grobner basis of /. Since y 2 x 2 — zx 3 sat-reduces to zero, G is the 
reduced Grobner basis of /. 

The number of S-pairs considered for reduction can be reduced drastically by using a version of 
the Gebauer-Moller criterion in algorithms (13.11) and (13.31) . The framework for properly explaining 
the Gebauer-Moller criterion is in the context of Grobner bases for modules (cf. [1J, §4). 
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3.1 The Gebauer-Moller criterion 

Let e\,...,e m denote the canonical basis of the finitely generated free module F = R m . A monomial 
in F is an element jc v e,, where x v is a monomial in R. Every element in F is a ^-linear combination 
of monomials. By definition a monomial x a ei divides a monomial x^ej if and only if i = j and x a 
divides x$ in R. We write this as x a et \ x^ej. A monomial order on F is a total order -< on monomials 
in F satisfying 

for every i,j = 1 , . . . , m and a,/3, y G N". We let in^ (/) denote the largest monomial in /. Now the 
Grobner basics for ideals in R can be generalized to submodules of F almost verbatim. For a subset 
B C F we let 'm^(B) denote the submodule generated by in^(/), where / G B. A Grobner basis 
of a submodule M C F is a set of elements G = {mi, . . . ,m t } C M satisfying in^(M) = in^(G). 
It is called minimal if in^ (m,) \ in^ (my) for i ^ j. We will use Grobner bases for submodules in 
reasoning about syzygies of monomial ideals. Consider a monomial ideal 

M= (x Vl ,...,x Vm ) CR. 

The syzygies of M are the relations in M i.e. the kernel K of the natural surjection R m — > M. Now 
consider the Z' 1 -grading deg(x,) = e ; on R. Then ^ is a homogeneous submodule of F in the Z7 1 - 
grading given by deg(e ; ) = v ; . A natural set of homogeneous generators are 

S tj = x Vi Vi v ■■(> ■, -x ViVv i- Vi ei G K 

for 1 < i < j <m (see \2\, Proposition 2.8). Define a monomial order -< (The Schreyer order) on F 
by x a ei -< x$ej if and only if 

a + Vi < P + Vj or a + Vi = /3 + Vj and i < j, 

where < is any term order on R. Then we have the following 

Proposition 3.5 The homogeneous generating set {Sjj \ 1 < i < j < m} is a Grobner basis for K 
over the Schreyer order -<. 

The Grobner basis in Proposition 13.51 is rarely minimal. In view of Theorem 2.9.9 in 0, it 
suffices to reduce the 5-pairs corresponding to a minimal Grobner basis of the syzygies (in Buch- 
berger's algorithm). This procedure is in fact one of the Gebauer-Moller criteria for cutting down on 
the number of 5-pairs. The point is that this minimization is easy and quite fast to perform in step 
(HieTl of Algorithm 13.11 Suppose that we must update S pairs with a non-zero polynomial g = g m , 
where G = {gi, . . . ,gm-l} m ste P <EH> - We put x Vi = m^(gi) for i = 1, . . . ,m. Consider the syzy- 
gies Sim, ■ ■ -,Sm-l,m- In the Schreyer order we have in^(5i m ) = * VlVv " , ~ Vm e m , . . . ,in x (S w _i )Hl ) = 
x v >"- lVVm ~ Vm e m . Thus the minimization can be done successively in step dneTl by throwing out super- 
fluous monomials among 

viVv m -v m 



^v m _i Vv m — v m 

This can be implemented as below (u < v means that v — uEN" for u, v G N"), where dnab represents 
the usual criterion, where leading terms are relatively prime (cf. Q, Proposition 2.9.4). 
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Algorithm 3.6 

updateSpairs: 

(i) MinSyz := 0; 

(ii) for each v ; - in {vi , . . . , v m _ i } do 

(a) if (v, n A Vi = 0) continue; 

(b) a = v,-Vv m -v m ; 

(c) if (w < a for some {w,p) G MinSyz) continue; 

(d) Delete {w,p) G MinSyz if a < w; 

(e) MinSyz : = MinSyz U { (a, S(g,-, g m )) }; 

(iii) for each (a,/?) G MinSyz do 

(a) Spairs := S pairs U {/?}; 

4 Lattice ideals 

Recall the decomposition of an integer vector v G Z" into v = v + — v _ , where v + , v~ G N" are vectors 
with disjoint support. For w,vG N" we let m < v denote the partial order given by v — u G N w . For a 
subset J C Z" we associate the ideal 

/^ = ( x v + _ x v- I v G ^) C 7?. 

In the case where ^ = jSf is a lattice we call 1% the lattice ideal associated to jSf , If w — v G Jz? for 
w,vG N", then 

x" -x v = x l '- {l '- v)+ (x {l '- v)+ -x (u - v) ~) G (1) 

The binomials B% = {x u — x v \ u — v G JSf } C are stable under the fundamental operations in 
Buchberger's algorithm: forming ^-polynomials and reducing modulo a subset of B _gf . This means 
that starting with a generating set for 1^ in 5^ we end up with a Grobner basis consisting of 
binomials in B _j? . Reducing a monomial x w by an element of B amounts to replacing x w by x w ~ v , 
where v G Jz? . Therefore if a binomial jc" — x v Eljf, then a — v G Jz? . This proves that / g» is saturated 
and algorithm (13.11) applies. The simple data structures in the specialization of algorithm (13.11) to 
lattice ideals are very appealing. If / = x" —x v , then 

sat(/)=^"- y ) + -x("- v ) _ . 

by (HJ). With this in mind we define 

bin(w) =x w+ -x w ~ 

for w G Z". Using this notation we have sat(bin(w),bin(v)) = bin(w — v). Similarly if v + < u + we 
may reduce bin(w) by bin(v). This results in a binomial / with sat(/) = bin(w — v). Notice that 



7 



replacing u by u — v if v + < u + corresponds to sat-reduction of bin(a) by bin(v). We have silently 
assumed that the initial term of bin(w) is x w for the term order in question. We will keep this 
convention throughout. 

Usually a generating set 3$ for Jzf as an abelian group is given. Computing the lattice ideal 
/jgf D lag can be done using that 

^jzf = Im- 

If SS contains a positive vector, then Ig = 1% (Q, Lemma 12.4). If J^f n N" = {0}, may be 
computed from 7^> using Grobner basis computations for different reverse lexicographic term orders 
(0, Lemma 12.1). 

With these conventions it is quite easy to convert algorithm (13.11) into a specialized algorithm for 
lattice ideals representing binomials via integer vectors with additional structure (like the degree of 
bin(v) and certain other (optimizing) features). We give the straightforward translation of algorithm 
(13.11) into the lattice case. 

Algorithm 4.1 (Homogeneous Buchberger algorithm for lattice ideals) 

INPUT: Term order -< . Integer vectors B = {vi , . . . , v r } with respect to -< such that (bin(vi ),..., bin(v r ) ) 
is a positively graded lattice ideal . 

OUTPUT: Integer vectors G = {w\,. . .,w s } such that (bin(wi), . . . ,bin(w,y)) is a minimal Grobner 
basis over -< for 1%. 

(i) Spairs:=<d;G:=<d; 

(ii) while (B ^ or Spairs ^ 0) do 

(a) Extract a binomial bin(v) of minimal degree in B U Spairs. 

(b) Compute the reduction bin(w) := bin(v) G ( sat ), continue if the degree drops in a sat- 
reduction step in the division algorithm. 

(c) if (bin(w) = 0) continue; 

(d) G:=GU{bin(w)}; 

(e) updateSpairs 



updateSpairs: 

(i) MinSyz := 0; 

(ii) for each bin(v) in G \ {bin(w)} do 

(a) if (w + A v + = 0) continue; 

(b) a = v+Vw + -w+; 

(c) if (u < a for some (u,p) E MinSyz) continue; 

(d) Delete (u,p) E MinSyz if a < u; 
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(e) MinSyz : = MinSyz U { (a, bin(w — v) } ; 
(iii) for each (a,bin(w)) G MinSyz do 
(a) S pairs := S pairs U {bin(w)}; 

Similarly algorithm |33] translates into 
Algorithm 4.2 

INPUT: Term order -<. Normalized integer vectors B = {vi,.. . ,v,} with respect to -<, such that 
(bin(vi), . . . ,bin(v r )) generates the ideal /. 

OUTPUT: Integer vectors G = {w\, . . -,w s } such that (bin(wi),. . .,bin(w 5 )) is a minimal Grobner 
basis over -< for / or proof that / is not a lattice ideal. 

(i) Spairs:=<d;G:=<d; 

(ii) while (B ^ or Spairs ^ 0) do 

(a) Extract a binomial bin(v) of minimal degree d in B U Spairs. 

(b) bin(w) :=bin(v) G ( sat ); 

(c) if (bin(w) = 0) continue; 

(d) if(deg(bin(w)) < d) 

(i) OUTPUT bin(w) as proof that / is not a lattice ideal and HALT 

(e) G:=GU{bin(w)}; 

(f) updateSpairs 

5 The Sullivant challenge 

Sullivant's challenge 4 is about deciding if the ideal J generated by a given set B of 145, 512 binomials 
generate the kernel P of the toric ring homomorphism 

k[x ijk ] ->k[u ih v ik ,Wj k ] 

given by xm i— > UijVi k Wj k , where 1 < i,j,k < 4. The 145,512 binomials are constructed by act- 
ing with a symmetry group on carefully selected binomials 5 . In this setting we need to compute 
in the polynomial ring k[xij k ] in 64 variables! The ideal J is homogeneous in the natural grading 
deg(xm) = ■ ■ • = deg(*444) = 1. The strategy is applying algorithm (14.21) to J using a reverse lex- 
icographic order. If algorithm (14.21) finishes without halting in step diidb . then Sullivant has proved 
that J must generate P. If not, algorithm (14.21) will halt with a binomial in P \ J. 

4 http : //math . berkeley . edu/~seths/ ccachallenge . html 

5 See http : //math . berkeley . edu/~seths/ccachallenge . ps for details 
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Running the gbasis command of GLATWALK with respect to the cost vector —e\ and the grading 

e\-\ h ^64 we compute a Grobner basis of J after converting the binomials in the two files 6 7 

containing J into integer vector format. After computing a 15-truncated Grobner basis, gbasis (in 
the incarnation of algorithm (14.21) ') outputs the degree 14 binomial 

2 

*3 1 1 *22 1 *43 1 *2 1 2* 1 22*342* 1 1 3*43 3*243*424* 1 34*334*444 ~ 

2 

*2 1 1 *42 1 *33 1 * 1 1 2*3 1 2*222*242*2 1 3* 1 33*443* 1 24*434*344 

as a binomial in J\ J proving that J does not generate P thereby answering Sullivant's computational 
commutative algebra challenge. Running gbasis in the above setting is not a simple computation. 
In fact the 15-truncated Grobner basis of J contains more than 300,000 binomials and the whole 
computation takes close to two days on most modern PCs. 

Details and more information, including the relevant files for Sullivant's challenge, are located 
at |http : //home . imf . au . dk/niels/ GLATWALK 
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